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SINGULAR  PERTURBATIONS  OF  ORDINARY  DIFFERENTIAL  E<5UATIONS 


Introduction 

Since  so  few  differential  equations  can  be  i)|tegrated  exactly, 
it  is  important  to  develop  approximate  methods  for  solving  such 
equations.  One  important  approximation  method,  and  almost  the  only 
one,  is  to  express  the  desired  solution  of  a  differential  equation  as 
a  modification  or  perturbation  of  the  solution  of  a  simpler,  that  is, 
more  nearly  solvable,  equation.  This  method  is  known  as  the  method 
of  perturbations  and  its  validity  usually  depends  on  the  fact  that 
the  modification  in  the  equation  or  in  the  solution  is  small  in  some 
sense.  Frequently,  the  smallness  of  the  perturbation  depends  on  the 
size  of  some  parameter  in  the  equation  and  v^e  may  be  able  to  express 
the  desired  solution  as  a  convergent  power  series  in  terms  of  the 
parameter . 

V/e  shall  discuss  the  perturbation  method  for  ordinary  differential 
equations  depending  upon  a  small  parameter  e .  We  begin  in  Chapter  I 
with  the  regular  case,  that  is,  the  case  where  naive  methods  give  the 
desired  solution  and  sometimes  even  give  it  as  a  convergent  power  series 
in  the  parameter.  As  illustrations  of  this  procedure,  we  shall  obtain 
the  Neumann  series  and  the  Fredholm  expansion  for  the  solution  of  .an 
integral  equation.  We  shall  also  indicate  how  to  obtain  the  solution 
of  the  eigenvalue  problem  for  an  ordinary  differential  equation  as 
a  power  series  in  e  . 


ii 


In  Chapter  II  we  begin  the  study  of  singular  perturbation 
problems,  that  is,  problems  in  which  the  standard  methods  fail  aind 
some  special  device  has  to  be  employed.  We  consider  differential 


■  .■3' 

equations  in  which  the  coefficient  of  the  highest  derivative  is 
multiplied  by  the  small  parameter  e;  therefore,  in  the  limit,  as  e 
approaches  zero,  the  order  of  the  differential  equation  is  lowered 
and  we  have  too  m^y  boundary  conditions.  The  question  of  whether 
the  limiting  problem  makes  sense  and -if  so,  which,  if  any,  of  the 
given  boundary  conditions  should  be  used  will  be  discussed  in  a 
very  general  case.  It  is  adso  shown  that  the  boundary  conditions 
which  are  lost  give  rise  to  a  "boundary  layer"  effect. 


In  Chapter  III  we  treat  the  relaxation  oscillations  for 
Van  der  Pol's  equation  with  a  large  value  of  the  parameter.  The 

problem  is  to  find  the  period  of  the  -  oscillation  in  terms  of  the 

# 

parameter.  ''Je  shall  obtain  two  terras  in  the  expansion  for  the 


period. 


Chapter  I 


REGULAR  PERTURBATIONS 

Iteration  and  the  Neumann  Series 

Consider  the  problem  of  solving  the  equation 
(1)  Lu-£Mu  =  f 

where  u  is  a  given  function  and  L  and  M  are  given  operators. 

We  assume  that  we  know  how  to  solve  the  equation 

Lu  =  f 

and  wish  to  use  the  solution  of  this  equation  to  solve  (1).  Suppose 
the  inverse  operator  L  ^  is  known.  Applying  it  to  (1),  we  get 

(2.)  u  -  eKu  =  g, 

where  we  have  put 

=  K,  L"^f  =  g. 

We  solve  (2)  by  an  iteration  method.  We  write  it  as 

(3)  u  =  g  +  s  Ku 

and  then  obtain  a  sequence  of  approximations  u^jU^jU^,...  by  putting 

(4)  u  =g+eKu,  n  =  0,l,2,... 

n+1  n  111 

Uq  =  g. 

There  are  now  two  problems.  First,  to  show  that  the  sequence  u^ 
converges  to  a  limit  and  secondly,  to  show  that  this  limit  satisfies  (3). 

The  usual  procedure  for  investigating  the  convergence  of  the 
sequence  begins  by  obtaining  an  estimate  for  the  difference  of 
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We  shall  not  discuss^  this  question  further. 

As  an  illustration  of  the  method,  we  consider  the  problem  of 
finding  a  solution  of  the  non-linear  ordinary  differential  equation 

(9)  y'  =  f(x,y) 

with  y(0)  =  y^.  Note  that  we  have  omitted  the  parameter  e  because 
it  is  not  needed  for  convergence  of  the  iteration  method.  Inverting 
(9) ,  we  get 

X 

y  =  Yq  +  ^  f[J^,y(?)]d?;. 

Just  as  in  (4),  we  set  up  an  iteration  procedure.  We  put  -  yQ 


and 

(10) 


'n+1 


(x)  =  y^  +  5  f[?,y  (i^)]d?^. 


For  a  norm  we  use 

Pn+1^^^  = 


'^n^l  -  ^n' 


ly,,l(x)  -y^(x)|. 


••  "  ■'  ^  sS 

-  tV  1»  X  *-< 

A  ■<  ^ 
A 


1 


To  prove  the  sequence  converges,  we  need  some  condition  on 

f(x,y).  We  assume  that 


9f(x,y) 


ay 


<  C 


for  all  y  and  0  <  x  <  a.  Then 


fC^,y^(^)]  -  f[^,y^_^(?)]  = 


(11) 


y  (O 


n 


9  T1 


dT] 


y^(0 


CdT, 


<  c|y^(!^)  -  =  Cp^(5), 


t 


.-‘'itn  * 


4 


.. 


".*4  *•  ,  4Vlf»>  A. 

■  ^  ■« 


;■---..'  ''  ■-. 


'i?'’ From '(10) ,  we  have 


using  (ll)'’we-get 

.  '.  •'  •  ■  „  •_  •  -x 


.■•/•:•■'••■  '•'•■VS; /(la)  '■••  <  c  ;  p„(?)dc. 


'•  ■  ■'•;•'■■  j  ;  ■■• 


;”  The  proof  of  convergence  can  be  completed  as  before  if  we  can 

obtain  an  estimate  for  p  (x)  such'  that  the  series 

,-.  •  n 

CD 

?  p  (x) 

:  ’1 

converges.  We  shall  show  later  how  such  an  estimate  can  be  derived 
from  a  more 'general  integral  inequality  than  (12). 


Let  us  consider  another  illustration  of  the  iteration  procedure. 

The  problem  is  to  find  a  solution  x(t)  of 

X  +  -x  .=  f(x,x,t) 

such  that  x(0)  =  A  and  x'(0)  =  0.  Here  we  may  consider  Lx  as  x  +  x 
ancf  Mx  =  f(x,x,t).  The  inverse  operator  to  L  is  easily  found  and  we 
obtain 

t 

X  =  Acost  +  i  f[x(T)  ,x(t)  ,T]sin(t  -  T)dT  . 

0  * 

We  set  up  the  customary  iteration  as  follows: 


Xq  =  A  cos  t 

t 

x^^^  =  Acost  +  j  f  [x^(t)  ,x^(t)  ,T3sin(t  -  T)dT  . 

» 

But  now,  we  see  that  the  estimate  for  x  , (t)  -  x  (t)  will  depend 
’  n+1  ^  n 

also  on  an  estimate  for  x(t)-x  -(t). 

n  n-1 


I 


5 


6 


where 

t 

M(t)  =  ]■  m(T)d‘t  . 

0 

This  suggests  putting  in  (14) 


p  (t)  =  «4-b  (t) 

n  nl  n 


After  this  substitution,  (14)  becomes 


(n  +  1)1  n+1 


7b  ,(t)  <  i  ^^^-b  (T)dT 

1  n+1  -  ^  nl  n 


Put  =  maix  b^(t)  for  0  <  t  <  a.  Then  this  inequality  implies  that 

(n  +  1)1  n+1  -n^  nl  n(n+l)l 


b  , (t)  <  B 
n+1  —  n 

« 

for  0  <  t  <  a:  consequently,  B  -,<B  <B  t<***<B^.  This  shows 

-  —  ’  n+1  -  n  -  n-1  —  —  1 

that  the  sequence  of  B  is  bounded  and  then 

n 

(15)  . 

Note  that  a  similar  result  can  b e  obtained  for  (12),  We  have 

m(t)  =  C,  M(t)  =  Ct  and 

^  „  (Ct)^ 

p  ( t )  ^  B  ,  • 

n  -  1  nl  s 

•  Using  the  estimate  given  in  (15),  we  find  that 

2  p  (t)  <  B  e^^^\ 

^  n  —  1  ’ 


a  convergent  result,  which  shows  that  the  iteration  procedure  defined 
in  (15)  converges  to  a  limit. 


The  iteration  method  leads  to  an  importaiit  result  if  the  operator 
l-(4)  is  linear.  In  such  a  case. 

? '.ViS  ■;  ;f  <■  •  ‘  - 


We  thus  obtain  the  formal  Neumann  serxes. 


The  series'  will  converge  .Vifi'the;  Operator, ;  .•:e.K  y.2is  .  .[.'small"  enough.  ■  A 


sufficient  condition  fqryicbnvergence  is  'that  the,  e.lgenvalue  of'  eK 


with  smallest  absolute  .value-. have:  magnitude,  less  than  one 


It  should  be  noted  that  ,  the.,  infinite-  series  (.16)  dould  also  be 


obtained  by  assuming  u  has  ..a  series  expansion  in  powers  -of  t.  such  as 


substituting  ■  this  in  (4.')  and  equating' the  coefficients  of  corresponding 


powers  of  e  on  each  s-ide  -of  the  '  resulting '  equation 


Fredholm  Expansion 


The  Neumann,  series  given  in 


sufficiently  small  values  of  e  .  -Fredholm  showed^hbw  to  obtain  an 


expansion  for  the  solution  of  (4)  which  is^'yalid  for  all  values  of  e 


for  which  the  solution  to  (4)  exists.  Of  course,  the  operator  K 


must  be  assumed  linear.  We  shall  obtain  this  expansion  in  a  formal 


manner  and  ignore  the  im  ortant  question  of  proving  the  convergence 
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The  solution  of  (4)  may  be  written  as  • 

(17)  u  =  (I  -  eK)~\ 

where  I  is  the  identity  operator.  If  we  consider  the  case  where  K 
is  a  finite-dimensional  matrix,  the  solution  (17)  breaks  down  when 
the  determinant  of  I  -  eK  is  zero.  Put 


D(e)  =  det(I  -  eK) . 

We  shall  state  later  the  meaning  of  the  determinant  of  a  non-matrix 
operator. 


(18) 

where 


Rewrite  (17)  as  follows 
eD(K,e)1 


u  = 


I  + 


D(e) 


g 


I  + 


eD(K,e) 
D(e  ) 


=  (I  -  eK) 


-1 


This  implies 


(19)  D(K,e)  = 


KD(e  ) 

I  -  eK  * 


If  K  is  a  finite-dimensional  matrix,  the  formula  (19)  suggests  that 
D(K,e )  would  be  regular  for  all  values  of  e  because  the  singularities 
of  (I  -  eK)  ^  occur  at  the  zeros  of  D(e )  and  would  be  cancelled 
by  the  factor  D(e)  in  the  numerator.  It  seems  reasonable  to  expect, 
therefore,  that  both  D(K,e)  and  D(e )  would  have  expansions  in 
powers  of  e  which  would  converge  for  all  values  of  e . 


Before  we  obtain  these  expansions,  we  wish  to  define  the  determinant 
of  a  general  operator.  If  e  is  small  enough,  we  may  define 
log(I  -  eK)  by  the  power  series 

-  ®  (ej^  , 

1 


o 


If  I  -  eK  is  a  matrix  with  eigenvalues  *  "  * ’^m’ 

lii 

det(l  -  eK)  =  X  \  =  exp  ( S’  logX  ). 

12’  ’  m  ^  y 

But  log(I'-  eK)  has  the  eigenvalues  logX^ , • • • , logX^  and 

m 

trace  log(I  -  eK)  =  E  logX  ; 

•  1  ^ 

therefore,  for  a  matrix  we  obtain  the  formula 


(20) 


det(I  -  eK)  =  exp[tr  log(I  -  eK)], 


where  tr  stands  for  trace. 

Formula  (20)  can  be  used  to  define  the  determinant  of  (I 
for  general  operators  K  if  the  trace  of  these  operators  can 
defined.  We  shall  defi^je.  trace  only  for  integral  operators. 


(21) 

then 

(22) 


Kg  =  5  k(x,y)g(y)dy 
a 


tr  K  =  i  k(y,y)dy. 
a 


Suppose 

CD  /  xn 

(23)  D(0  =  2  Vr- = 

„  nl  n 


and 

(24) 


CD  ,  ,.n 
D(K,e)  =  s 

O  n.  n 


where  c^  are  constants  and  are  operators.  Using  (20), 


that 

(25) 


CD  ,  .n  CD  n 

log  D(e)  =  trClogd  -  eK)]  =  -tr  E  iiiii-  =  _  ^ 

^  n  ^  n  n 


:  -  eK) 
be 

Suppose 


we  find 
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where 


k  =  tr  K 
n 


Differentiating  (25),  we  get 


2  e^-^k  -  -  tr  ^ 

D(  e )  •  ^  .  n  "  ^  I  -  eK  ’ 


OD  .  vn  ®  1 

D'(e)  =  2  \  ,  c  £  =  -  (r  e  k  )D(e), 

^Cn-Dln  ^  n 


Multiplying  the  two  series  on  the  right-hand  side  and  equating  coefficients 

of  corresponding  powers  of  e,  we  get  a  set  of  linear  equations  from 

which  c  could  be  determined.  Soon,  we  shall  obtain  k  more  useful  • 
n  ’ 

formula  for  c  . 

n 

From  (19)  we  have 

D(K,e)  -  eKD(K,e)  =  KD(e ) . 

Substituting  the  series  (25)  and  (24)  in  this  equation  and  comparing 
the  coefficients  of  corresponding  powers  of  e ,  we  get 

(28)  B  +  nKB  ^  =  c  K 

n  n-1  n 

for  n  =  0,1,2,....  Note  that  from  (25),  c^  =  1,  =  K. 

Again  use  (19)  and  taike  the  trace  of  both  sides.  We  find 

tr  D(K,e)  =  D(e  )  tr  ^  „  =  -  D' (e  ) 

by  the  use  of  (27).  Comparing  corresponding  powers  of  e ,  we  obtain 
the  promised  formula  for  c^,  namely, 

(29)  c^  =  tr 


I 
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Formulas  (28)  and  (29)  enable  us  to  determine  c  and  B  by 

n,  n 


succession.  We  find 


(30) 


<=1  =  ,  \  =  k^K  -  K 


'^2  "  ^1  ■  ^2’  ®2  "  ^^1  "  ^2'^^  ~  ^^®1’ 


etc.  This  is  as  far  as  we  can  go  in  the  general  case.  However,  if 

O 

we  assume  that-  K  is  an  integral  operator,  as  Fredholm  did,  we  can 
.give  some ..elegant  determinaptal -representations  for  and  c^. 


We  .assume 


.  ■  -.Kg  =  j,  K(x,y)g(y.).dy 

’  .  .  a  .  '.  '  . 

where  K(x,y)  '  is  called  the  kernel-  of  the-  integral  operator  K. 

■  If  B  (x,y)  is  the  kernel  for  B  ,  ■  then  (28)  becomes 
-  .  .  n-  ’  .  •  .  •  n’ 


(31) 

and  .(29)  becomes 


Bn(x,y)  +  n  ^  k(x,z)B^j_^(z,y.)dz  =  c^K(x.,y), 
'  '  a  . 


(32)  =  5.  B^_^(y,y)dy  . 

a 


With  these  notations,  the  formulas  (30)  may  be  written  as  follows: 


b  . 


=  *  Bt'l 

a 


K(to,t^) 


K(t^,t2) 


K(t^,t^) 


=1=1  K(t^,t^)dt, 

a  "  • 


•  1 


m 


12 


b 

a 

It  is  clear  what  the  general  form  of  B  (t  ,t  )  and  c  should 

“  n  0  n+1  •  n 

be.  By  mathematical  induction  using  (jl)  and  (32),  it  is  not  difficult 
to  prove  the  formulas  correct.  We  leave  the  details  to  the  reader. 


b  b 

I  5  dt^dt^ 


KCtQ.t^)  K(t^,t^) 
K(t^,t^)  ®K(t^,t^) 


KCt^.t  ) 


K(t2,t^) 


K(to.t2) 

K(t^,t2) 

K(t2,t2) 


Perturbation  of  Eigenvalues 


Consider  the  problem  of  finding  the  eigenvalues  of  the  linear 
operator  L  +  eM  if  the  eigenvalues  and  eigenfunctions  of  L  are 
known.  The  problem  is  to  find  values  of  X  ■  for  which  there  exist 
non-zero  solutions  of 

(33.)  (L  +  eM)u  =  Xu. 

Here  both  L  and  M  are  linear  operators  ancf  we  assume  L  is 
seif-adjoint. 


Suppose  Uq 
eigenvalue  X^ 
approaches  X^ 


is  an  eigenfunction  of  L  corresponding  'to 
and  we  wish  to  find  the  eigenvalue  of  L  +  eM 
as  e  approaches  zero .  Assume ' 


the 

which 


and 

CD  k  •  '  '  ’ 

u  =  u_  +  2  e  u, 

0  ^  k 

where  the  X  and  u  (k  >  1)  are  to  be  determined.  Substitute  these 
k  k  — 

expansions  for  X  and  u  into  (33)  and  equate  the  coefficients  of  the 


corresponding  powers  of  e  on  each  side  of  the  resulting  equation. 
In  .this  way  we  obtain  the  following  set  of  equations: 


(34) 

%  =  Vo- 

(35) 

(L  -  4q)u^  = 

(Xi 

-  M)Uq  • 

(36) 

etc . 

(L  -  ?^q)u2  = 

(Xi 

-  M)u^  +  Vo’ 

® 

The  first  of  these  equations  is  automatically  satisfied  by  our 
choice  of  and  u^.  In  the  second  equation  we  do  not  know  X^ 

and  u^ .  In  general,  this  equation  has  no  solution  because  the 
homogeneous  equation 

(L  -  Xq)v  =  0 

has  a  non-zero  solution  v  =  u^.  A  necessary  condition  for  the 
existence  of  a  solution  of  the  second  equation  is  that  the-. right- 
hand  side  be -orthogonal  to  that  is, 

i'uo(X^  -  M)Uq  =  0.  . 

This  equation  implies 
•  .  ^1  ,-^^0  =  ^ 

and  thus  defines  X^ .  With  this  definition  of  X^,  equation  (35) 
will  in  general  have  an  infinite  number  of  solutions  because  to  any 
particular  solution  of  (35)  we  may  add  where  is  an 

arbitrary  constant.  We  choose  so  that  if 


then 
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Chapter  II 
/ 

SINGULAR  PERTURBATIONS 


Introduction  *  ■  .  .  ■ 

In  the  last  section  of  Chapter  t  we  have  discussed  a  method  for 
finding  the  eigenvalues  of  a  complicated  operator  as  perturbations  of  the 
eigenvalues  for  a  simpler  operator.  In  the  present  chapter  we  shall 
discuss  operators  for  which  .the  proposed  method  breaks  down  and  for 
which  a  more  sophisticated  approach  Is  needed.  The  type-  of . difficulty  ■ 
we  shall  treat  is  Illustrated  in  the  following-  example  given  by 
Lord  .Rayleigh: 

Find  the  values  of  \  for  which  there  exist  non-zero  solutions  of 
(1)  .  -  .  eu^^  -  u"  =  \u‘ 

such  that  u(0)  =  u'(0)  =  u(l)  =  u’(l)  =’0. 


This  eigenvalue  problem  when  treated  by  the  method  of  the  preceding 
chapter  leads  to  the  equation  .  •  ,  -  '  • 

(2)  •  ■  -  u"  =  Ku 


with  the  boundary  conditions  u(0)  =  u'(0)  =  u(l)  =  u'(l)  =  0.  For  any 
value  of  the  only  solution  of  (2)  satisfying  all  four  boundary 

conditions  is  u  Identically  zero  and  the  perturbation  method  breaks 
down  before  we  can  even  begin. 

A  problem  such  as  (l)  for  which  the  simple  perturbation  techniques 
explained  in  Chapter  I  do  not  work  will  be  called  a  singular  perturbation 
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problem.  The  singularity  is  obviously  due  to  the-  fact  that  the  zeroth 
order  equation,  that  is,  the  equation  for  e  =  0,  is  of  lower  order 
than  the  original  equation  (l)  and  we  are  troubled  by  a  plethora  of 
boundary  conditions.  This  chapter  will’ treat  the  general  case  of  an 
operator  L  +  eM,  where  L  is  a  differential  operator  of  the  nth 

order,  M  a  differential  operator  of  the  mth  order  and  m  >  n.  We 

* 

shall  discuss  both  the  eigenvalue  problem  for  L  +  eM  and  the  solution 
of  the  non-homogeneous  equation 

(L  +  eM)u  =  f. 


and  also  give  a  priori  niles  for  determining  which  boundary  conditions' 
will  be  used. 


■  The  Exact  Solution  for  the  Rayleigh  Problem  .  ' 

•  Because  equation  (l)  has -.constant- coefficients ,  it  can  be  solved 

'exactly..  Assiime.  u  =  e^^"  is  a  s'olu’tioh  of  (l) }  -then  r . ,  must'-satisfy 

the  equation  ‘  . 

-  ■■  ■  ■  ■  •  4  ’  2  ■■■■ 

•  ■  er  -  r  =  X. 


Solving  for-  r  ,  we  get 
(3) 


2 

r  = 


1  +  Vl  +  4Xe: 
~2e 


We  assume  e  '  small  and  positive  and  equal  to  n  and’conslder  first  the 
.positive  square  root-  in  (3) .  We  find 

^2  2.+  2Xny.,pi.ri^)  ^  ^-2^^  ^  ^^2  ^  j  ^ 

2n'' 
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the  limit  we  must  have 

sin  =  0. 

However,  from  (5)  ,  ^2  ~  limit;  therefore, 

(7)  sin  J\  =  0, 
or  K  =  n^TT^  for  n  >  1. 

•  (S 

This  result  shows  that  in  the  limit  the’  eigenvalue  problem  for  (l) 
goes,  into  the  eigenvalue  problem  for  (2)  with  the  boundary  conditions 
u(o)-  =  u(l)  ■=  0.  Using  (7)  and  the  boundary  conditions,  we  can  see  that 
also  the  eigenfunction  of  (l)  approaches  the  eigenfunction  of  (2)  with  '  ’  ■ 

zero  boundary  conditions,  namely,  sin  nnx. 

Boundary  Layer  Effect  •  •  °  ' 

We  have 'just  shown  that  as  '.e  .^.pproaches  zero  the  eigenfunctions 
for  (1)  approach  sin  nrrx.  ’  This  seems -paradoxical  -because,-  even  though 
u  =  sin  niTX  satisfies  the  lirai'ting  equ’ation.  (2)  and  some  of  the  boundary  • 
conditions,  it  does  not  satisfy  the ’  conditions  -’u’CO)  =  u'Cl)  =  0--  •  We  . . 
proceed  to  investigate  how  the  boundary  condition  u'(0)  0  is  lost  '  .  •-  • 

as’  s  approaches  zero.  •  ’-•  • 

To  study,  the  behavior  of  the  solutions  -of  (l)  near  ’  x  =  0,  ’we 
"stretch"  the  x-axis  by- putting  -  •  '  '  ^ 

X  =  Tit  '  ■  •  . 

in  (1) .  The  equation  becomes 

(8)  -  u  =  \q^u 
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where  dots  denote  differentiation  with  respect  to  t.  The  boundary 
conditions  at  t  =  0  are  still  u(0)  =  u'(0)  =  0  but  the  boundary  con¬ 
ditions  at  X  =  1  become  conditions  at  t  =  f]  ^  which  approaches  infinity. 
Notice,  however,  that  t  =  oo  does  not  correspond  to  one  value  of  x  but 
to  all  values  of  x  such  that  f)  ^x  is  unbounded.  Because  of  this,  we 
cannot  apply  the  conditions  at  x  =  1.  Instead,  we  try  to  find  a 

solution  of  (8)  or,  more  precisely,  of  the  limiting  form  of  (8) ,  that  is, 

(9)  u*-u  =  0, 

such  that  u  will,  for  large  values  of  t,  fit  in  or  match  the  function 
Uq  =  sin  rnx  =  sin  TrnTjt.  ^ 

Any  solution  of  (9)  that  satisfies  the  boundary  conditions 
u(0)  =  u'(0)  =  0  can  be  written  as  .  ■ 

(10)  u  =  a(l  +  t  -  e^)  +  p(l  -  t  -  e"'^)  -  .■ 

where  a  and  p  are  arbitrary  constants.  To  have  this  ’solution  fit  into  • 
the'  function  sin  nnrit,  it  is  clear  that,  first,  a  must  be  zero.  ■  Next, 
if  we  let  ■  t  ■  go  to  infinity  but 'in  such  a  fashion  that  x  is  small  (this, 
can.be  done  by  taking  ■  t  =  r]  '  .  'and  x  =  B  '  ■') »  "the  function  /  sin  nniit 
can  be  approximated  by,  TTnT]t;'  consequently,  for  u  to  behave  like 
.irnrit  when  t  is  large,  we  must 'take  p  =  --TmT].  We  therefore  conclude, 
that  the  appropriate  solution  of  (9)  is  ■  .  ,  • 

(11)  u  .=  TTnqCe  ^  -  1'+  t)'  ■=  nnx  +  TTnTi(e  -  l)..  •' 

A  similar  argument  can  be  made  for  the  neighborhood  of  x  =  1.  If 
we  introduce  the  stretched  coordinate 

1  -  X  =  TIT 


into  (l)  and  go  to  the  limit  f)  =  0,  we  again  obtain  the  equation 

(12)  u^^  _  u"  =  0 

with  the  conditions  u  =  u'  =0  for  x  =  0  and  u  behaves  like  sin  nnCl-qx) 
as  X  approaches  infinity.  The  solution  of  (12)  satisfying  these  con¬ 
ditions  is 

(13)  u  =  cos  nn  (1  -  x  -  e~'’')T7nri 

=  cos  nTT[TTnf)(l  -  -  TTn(l' -  x)]. 

A  study  of  (ll)  and  (13)  explains  the  paradoxical  behavior  of  the 
zeroth-order  elgerifunctions  at  the-boundary  conditions.  ■'  The  actual  eigen-' 
function,  for  t)  0,  behaves  like  (ll)  in  the  neighborhood  of  x  =  0 
and  like  (13)  "in  the'  neighborhood  of  x  =’l;  thus the  actual  eigenfunction 
■satisfies  .all  four  boundary  conditions.  'But  as  both  (ll)  and  (13)  indicate, 
as  'x  moves  away- from  either  endpoint’ the  exponential  terms  damp  out  very 
quickly  and  the  actual  e'igenfViniction  goes  .over ’into  sin  nnx,  the  .zeroth- 
order  eigenfunction.'  ■  This  ’ situation  is  illustrated,' -on  an  exaggerated 
scale,  in  'Flg.  1.,  ^  ■  ■  ’■  ..  ..■  .  ■  ■ 


u 


FIGURE  1 


The  behavior  of  the  curve  at  x  =  0  and  x  =  1  is  typical  of  that  which 
occurs  in  boundary  layer  problems  in  fluid  flow.  We  shall  call  bovmdary 
layer  effect  any  such  occurrence  in  which  a  very  sudden ' transition  occurs 
between  two  parts  of  a  solution. 
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First  Order  Correction  to  the  Eigenvalues 

To  obtain  an  improved  value  for  the  eigenvalues  of  (l) ,  we  use  (8) 
after  replacing  u  on  the  right-hand  side  by  the  zeroth  order  approxi¬ 
mation  sin  X  X  thus  obtaining 

(14)  *u*  -  u  =  f]^Xsln^/X  T)t  . 

We  solve  this  equation  in  the  neighborhood  of  x  =  0  with  the  boundary 
conditions  u=u'=0  at  t  =  0  and  u  matches  sin  nnTit  as  t 
approaches  infinity.  A  particular  solution,  of  this  equation  is  the 
function  _  .  ,  .  . 

■  sinj\  pt 
2 

.  1  +  fl'^X 

•For  the  complementary  solution  we  shall  use  the  fractions  e"*^,  e  and 
the  functions  cos  liTrqt'  and  sin  nTiqt  instead  of  the  functions  1  and 
t.  We  choose  the  cosine  and  sine  functions  because  they  are  the  exact 
solutions  of  the  rimlt  equation  (2).  We  therefore  put 

u  =  sin  ^/X  T]t  ^  nTTT)t  +  p  sinn-rrqt  +  Y_®  ^  ' 

1  +  ri  X  .  ; 

where  a,-  p,  y  ®^®  't®  t>e  determined  to 'fit  the  boundary  conditions. 

In  order  that  u  not  go- -to  infinity  as  t  goes  to  infinity,  we  must 
have'  “  Q*  To  satisfy  the  conditions  u  =  0  and  u'  =0',  we  find  • 

•  .  ■  a  +  Y_  = 

n/x  T)(P  +  - -  Y  =  0. 

1  +  T]  X 
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To  have  u  match  with  sin  nirnt  as  t  goes  to  infinity,  we  must  have 

p  +  (1  +  =  1, 

and  spK  approach  nrr.  Neglecting  terms  of  order  we  conclude  that 


(15) 


u  = 


sin 


-  «yr  Ti(cos  nTTTit  -  e~^) 


1  + 

sin  jc 

1  +  f,\ 


-  Tl(cos  mrx  -  e  . 


Formula  (15)  gives  the  correct  (to  order  f]  )  behavior  of  the  eigen¬ 
function  near  x  =  0  and  also  in  the  middle  of  the  interval  0  <  x  <  1. 

To  find  the  correct  behavior  of  the  eigenfunction  near  x  =  1,  we  put 
1  -  X  =  r]T;  in  (l)  and  obtain 

••••  —  . 

u  -  u,.=  \U.  ■ 

We_de_sire  the  solution  of- this  which  satisfies  the  conditions  u  =  u'  =  0  . 
at  T  =  0  and  again  matches  sin  n/\  x  =  sin  (  1  -  tit)  as  1-  approaches 
infinity.  Just  as  in  '  (15) , .'we  . find  that  ' 

—  sin  n/X  nx 
■  .  ■  ■  cos  n/\ 

(16)  •  .  .  .  ■ 


■  .  •  ■  cos  'J\  (1+ri^)  .  , 


+'  r)  (cos  nirriT  -  e~^)’ 


sin  ►/X  (1-x)  ,  ■  /r  r  \  ■  -(l-x)/m 

=  - ^ ^ +.  Ti  n/.X  [cos  njT(l-x)  -  e-  '  "  IJ 


(l+Ti  X)cos.VX' 


correct 'to  terms  of  order  t] 


■ In  order  that  (15)  -and  (l6) '  represent- the  same  solution  we  must  have 


sin 


X  —  Tl 


T)  is/\  cos  niTX 


or 


(17) 


_  sin  n/X  (1-x)  ,  /r 
=  .  - — <-  +  Ti  n/X  cos  nr  cos  niTX, 

cos  n/X 


cos  is/x  X  sin  iv/X  =  t]  n/X  cos  nTTx(l  +  cos  cos  nir)  • 
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power' series  in  t]..  (it  Is  assumed  that  P^(x)  and-  q^(x)  are  analytic 

functions  of  x  in  a  <  x  <  b  and  0  in. that  ’interval.)  '  We 

indicate  the  idea.  Let  -u^^Cx)  ,be  any  solution  of 

Lu  —  0 

neglecting  boundary  conditions.  Put 


.(x)  =  u_(x)  +  2  Ti  u.,(x), 


1  <  j  <  n 


in  (21)  arid  compare  corresponding  powers  of  t).  We  find 


Lu.j^(x)  =  -Mu.^j^_^(x). 

Let  u.,(x)  be' the  solution  of  this  equation  such  that  the  solution  and 

JK  .  ,  •  .  ■  • 

its  first  n-1  derivatives  are  zero  at  x  ==  a.  It  is  clear  that  this 

sol-ution  exists  and  is  unique.  Thus,'  starting  with  u._(x),''  we  can 

J  ^ 

'successively  find  all  the  u  (x)-  •  and  thus  obtain  the  expansion  (22)  for' 
.  .  -  ■  ■  ■’ 
'u.(x).  It  can  be  shown  that  the .expansion  is. an  asymptotic  one.  ^ 

choosing  ■u^q(x),  .(1  <  j  <  n)  as  n  linearly  independent  solutions 

■of  .Liu'  =  0,.  it  is  clear  'that  we-  shall  get'  n  linearly  independent  • 

solutions .  of  (2l)  ..  ■  ' 

.  ,  To  obtain  t  other  linearly 'independent  solutions  of  (21),  we  begin' 

by  considering  the . characteristic  equation  ' 

(23)  Pj^(x)  +  r\^(x)  =0. 


If  we  assume  f'  0  <  x  <  b,  this  equation  has  for  every  x 

in  (a,b)  t  solutions,  namely,  the  t  t  ^  roots  of  -p^(x)/q^(x) .  We 
shall  assume  that 


Re  r(x)  f  0 


for  every  x  in  (a,b)  and  any  root  of  (23) 


For  any  root  r(x)  of  (23)  we  shall  construct  an  asymptotic 
expansion  of  a  solution  to  (21).  For  simplicity  of  exposition,  we  shall 


show  how  the  first  two  terms  of  the  expansion  can  be  obtained.  The 


infinite  expansion  can  be  obtained  in  a  similar  fashion.  For  further 


details  we  refer  to  a  paper  by  Jurgen  Moser  in  Communications  on  Pure  and 


Applied-  Mathematics,  vol.  VIII,  1955 i  ;pp.  251-279 


It  is  important  to  distinguish  between  Re 'r(x)  >  0  and  Re  r(x)  0 
We  shall  denote  a  solution  of  (21)  corresponding  to  -the  case  Re  r(x)  >0' 
by  .y(x) :  and  a  solution  corresponding  to ''the  case  Re  r(x)  <0  by  ■w(x). 


'h(,x).;;  =  '  5  r(  O.d^ 


where  '  :.a. 


Note  that 


This  implies  that 


o 
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.  -  n  ,  m\  ,  ’0^ '  r  /  -1  %  n-2  ,  ,  ,  v  ni-2T  ,  r  '  n-1  , 

■Vi.(p^r  +  qjj^r  )  +  — ^[n(n-l)p^r  +  m(m-l)q^r  ]  +  + 

,  ■  in— 1 T  .  ,  ,  /  n  —1  ,  m— 1  \  

•■"W  )-0. 

Because  of  equation  (23),  the  first  of  these  equations  is  satisfied  for 
arbitrary  ^^(x)  and  the  second  equation  reduces  to 


(24)  Vq  •  [np^(  x)  r’^  ^  +  mq^(x)r'^~^]  +  Vq[  2 


n(n-l)  n-2  ,  , 
— ^ or  r'  + 
'■^n 


••  + 


m(m-l)  m-2  ,  .  n-1  m-l-i  _  „ 

— — -  qr  r'+p,r  H-q^r  J-0. 
2  ..  m  •  •  -^n-l  ^m-1  ^ 


The  coefficient  of.  Vq'  in  (24)  can  never  be  zero  because,  if- it  were, 
r  '  would  be  a  multiple'  root  of  (23)  —  which  is  known  'to  be  impossible; ' 
consequently,  (24)  is  'a.  linear  ordinary  differential  equation''for 
which  can  be  solved  by ■ separation  of  variables. 


.  It  .is  easy  to  show  that  a  similar  procedure- will  enable  us  to  I,  -..- 
■determine. a  complete  expansion  for  v(x) .  The -proof  that  it  is  an  '  ' 
asymptotic  expansion  of  a  solution  of  (2l)  will  be  found’'in  the,' paper  of 
J.  'Moser  cited  previously..  )  '  .  _ 


The  Lost  Boundary  Conditions  ...  ,  ’  •  ' 

Let  us  return  to  (18) .  For  the  non-homogeneous  problem  let  .u^  be 
a  particular  solution  of-(l8).  Then  a  solution  u(x)  of (21)  must  be 
found  such  that  u+Uq  satisfies  the  boundary  conditions  (19) .  Similarly,  " 
in  the  eigenvalue  problem  a  solution  u  of  (21)  must  be  fo'und  satisfying 
the  boundary  conditions  (19).  We  wish  to  determine  which  boundary  conditions 
will  be  retained  in  the  limit  e  or  t]  equals  zero. 
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The  previous  sections  have  sho-wn  that  the  solutions  of  (21)  can  be 
divided  among  three  classes:  the  functions  Uj^(x) , . . .  ,u^(x)  which 
were  obtained  in  (22)  by  starting  ^d.tK  solutions  of  Lu  =  .0,  the 
functions  v^(x) , . . .  ,Vp('x)  ,  where  p  is  an  integer  not  greater  than  t, 
which  were  obtained  in  the  preceding  section  on  the  assumption_  -Re  r(x)  >  0, 
and,  finally,  the  functions  w, (x),...,w.  (x) ,  which  are  obtained 

X  -  Li 

similarly  to  the  v(x)  but  on  the  assumption  that  Re  r(x)  <  0. 


If-  u  is  a  solution  of  (21) ,  we  may  write 


u  =  u^  +  ••*.+  a  u  +  6^  V-,  +  •  •  •  +  B  V  +  w^  +  •  •  •  +  Y ,  w , 

11  nn^^ll  '^pp'll  ‘  t-p  t-p 

arid  then  the 'boundary  conditions  will  .give  us  equations  such  as  these: 


(25)  ■  ■  ■  •  -  ■ . 


m. 


The  existence  of  a  solution  of  these  linear  eguatlong  will  depend  bn  the. 
behavior  of  the  determinant  of  the' coefficients,  that  is',  -on  'the  determinant 


V^i) 

B^(uJ 

.Bi(vi),---. 

.Kl(Wi) 

■ • •  . 

.^..■B^(Ur^): 

■■■ 

.■■Bs(wi) 

•••  y”t-p> 

®s+r^^i^ 

'•  •  •  B  ,  (u  ) 
s+1  n 

.Vi^y 

®s+'l^^l^ 

...  B  , ,  (w,, 
s+1'  t-p 

•  • •  B  (u  )  . 
m  n 

Bjvi)  ... 

• ••  B  (w.  ) 

m'  t-p 

Suppose  that  ■"  v. 


where 


and  Re  >0;  then  .the  quantity  exp[t^’’  n,  (b)  ]  will  be  expbnentially 


If  we  multiply  ,  the  column  pf  .■■-,A,.,'^’'cbh' 


the  first  s  row.s  will.;  have';.- the;:' factor 


mm 


which  goes  exponentially  .to  .zero  .as  f)  '.  goes  to'.pero;-;  Thi's '.implies '.f hat 


the  first,  s  rows,  of  the  v-columns  will  go  to  zero 


Similarly,  by  multiplying  the ..  w-p plumn's  ,  _-by 


last  m-s  raws', 'will,  have  .‘the  factof 


exp  [  ti“  h(b)  h(a)  ]  exp['-ri";:...  fdcr(- 


which  goes  exponentially  to  zero  as 


;This  impli‘es;.:.,thaf 


the 'last'  m-s  .rows  of '  the'  w-columns.. 'will  gp/tbv.zefo 


There,  is  one  .further  fact..'-  ..to  ■..'be  ;'n'dticed 


of  the  exponential.  As.  r|  goes-'to  zerc)",  ■  this  ■  implie's''..that  ’''.  B 


goes  to  zero  compared  to. 


to  zero  compared  to  '.for  '.1  C  k„' <;m-s-l'i  'Using  ,these'\;fac'ts;.,wev;' 

see  that  the  .leading  term  in,  ''A  is  obtained  by  expanding  the '■.folio'wing- 


determinant 


at  the  other  endpoint] 


«  Discard  the  boundary  conditions  in  decreasini 


order  of  magnitude  of  the  highest  derivative  present  in  the  boundar; 


condition. 


.(26) 


••• 


0  •••  B^«(w  ) 


. v.:'dierei;the ’'primed:!  <;B -i'' -an^',  B-.  i  and  B  modified  by  the  - 

.'t '  ■«■■■•■' elimlnatioH' of  the  •  exponential  factor  and  the  power  of,  t].'  Clearly,  the 

this'.'de'ternilh'ant  will  -not  contain  any,  term  involying  ■ 


for  -  the' 


'•t  ’ ' will -not  contain  any ,  term  involying  ■  ' 

t-  ('^v)  ■•'.'  ■ '..Tikis’  means ■''the''''bO'^^  -  Of  .highes.t  derivative  -• 

-V..J  ' ' '  ^"'1 "  '  '  '  ■■•.  ■'  ■V.v,"'  *■  -  “  *i4 i,'  '  ■  .'■  . ■’*'■.  ■  -  ■ 

.»  -  •  '  !,  J'-  ■■  .  y-'  ■  .  ‘.'s.'...  '  '  V'.‘  ■ 

.dS.',  •  ;  '.border';  atli, a.,  ■'and  ,V  b  V  niust-.be'.'ciiscarde^^  only  ■' m-2  conditions  are  t'etained.' 

"ttTh4,-4®'terM^nant';iht't^  zero '.because'.'there.  are  ,,.  '  ,' 

^'''-'fin-tioo  mahylzero.  rows'V7"  Thisriwili.'always'-hap.pen'/l-f7-.k'<-.  mj-.2..  .'in ’-that'  case.  '  ’  . 

1  1:,.  .  ■■.'•■■■•'.f--','' .A,';.-,  'f.'' ih.'the  ■  expansion  .'of '''A.,'' we''"inust  retain  row  '.ni-l.and  tow  ■•m-1  for -the' 


'  .  '  :  .  .Y'sI  ,.and.J^rpwi  s 'land,;, .row  1;  st-l"'^  ..the'-^'.'w'.s;  Thislwould  imply  .discarding  .  ,;  ’■ 

iViif.-:'  ■j.>'‘'Yithe''’tvd'.highes'£ 'deriyati.ve ’b6undaiy'''cpnditions  .■■for- bpjh  '  x'=.'a  -.  and  x.=  b. 

■  This?  procedure ^contlnuesfuntil  .we? obtain  a  ,nbh-idfehtically  zero  .determinant  ■ 

wlll'-'be  equivalent, ,  exactly  n;  boundary  cpnditlohs  -retained. 

■  '-d-  '  ’  ' -An -pxamlhation  of '  this'-’raethdd  leads  to  .the,  following  rule;  -  , 

■■'■  Discard  a  boundary- condition  at  an  endpoint  for' each  solution 


of  (21)  that- becomes  eOTonentially  large  there  (relative  to  its . behavior 
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For  the  present  problem  since  -there. .are  p  functions  v(x)  which 

become  exponentially  large  at  x  =’b,  we  discard  boundary  conditions 

B  B  and  since  there  are  t  p  functions  w(x)  which 

m  m-1  '  m-p+l  •  ^  ' 

become  exponentially  large  at  x  =  a,  we  discard  boundary  conditions  B  , 
B^  •*•»  t+p  1*  count  is  right  because  t  conditions  will  be 

discarded  and  m  --t  =  n  will  be  retained. 

■  What  happens  if  there  are  more  exponentially  large  functions  at  an 
endpoint  than: there  are  boundary  conditions?  For  example,  what  happens 
if  s  -  t  -f  p.<  1?  In -that  case,  the  solution  of  (18)  will  not  approach 
a  limit -as  e  ■  gpes  to  zero.  Consider  -the  following  example,:  ■  .  ' 


.  (27)- 


’eu"  au'  =1 


with  uC'O)  =.  u.'(O)  =  0.'  The  solution  of  this  equation  .is 


^  -  a--'  -:? 
■  • •  a 


(-1  + 


If a- >  0,  so  that  the  'solution  e  .'  is  exponentially  .large  at-  x'-= '0 
■(compared  to' its  valRie  e  ’’at  .pc  •=  l)j.  then  (28).  approaches  the  limit 
..  x/a.  This;'- limit  is  also  the  solution. the- limit  of  equation  (27),  . 


namely,  ■ 


au '  —  1 


with' the- highest  deri-vative  boundary  co'ndition  (u' 0)  discarded  and 
the'  other  boundary  condition  -  (u  =  O)  '  retained.  ",  .  , 

'  ■  '  ■  '  ■  ■  'i'ax/e  ''  'c'  -  - '  "  ■ 

If  a  <  0,  the  "solution,  e  '  ". -is"- exponentially-, large  at  x  =  1;  ' 
but  there  is  no -boundary  condition  at  x-  =1.  Note,  .ho.wever,  that, (28) 
doe.s  not  converge  to  a  limit,  in.  this  case. 
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.  •  ■  .  .  Chapter  III  ^  .  .  ••  ••  ;  .  ■  .  ; 

•  ■  '•  RELAMTION  OSCILLATIONS  .  '  '  '  .  . 

Ati- interesting  type  of  singular  perturbation  occurs  in  the  study  of 
the  Van  der  Pol  equation 

(1)  y”  -  p(l  -  y^)y'  +  y  =  o 

® when  the  parameter  ^  is- very  large.  If  we  change  the  independent’ 

—2  • 

variable  from  t  to  pt  and  put  e  =  p~  ,  equation  (l)  becomes 

(2)  ey  -  (1  -  y2)y  +  y  -  0 

.  \9  .  #1  ■  ' 

where  dote  denote  differentiation  with  respe’ct  to  the  new  variable.  :  Put 
y  =  x-  in  (2)  and- integrate  to  get  the- following  form  of  the  equa-tion: 

.  -  ■  ■  -.5  .  .  ■  >  - 

.  (3)  •  .  ,€X  -  (x  -  ^)  +  X  =  0. 

■  ■  We  shall  show  that  (3)' has  periodic  solutions  and  we  shall  find' the 
period  T  asymptotically  in'  e  for  small  values  of  e.  From  the -form’  - 
of  (3)  it  is  clear  that  for-  e  =  0  the  differential  equation- reduces  to 
a  first-  order  equation  which  has  no  periodic,  solution;:  consequently,  the 
.’perturbation  around  e  =  0  is  singular.-  ■  ’  .  .  -  - 

The  analysis  of  (3)  becom'es  simpler  when  time  is- .eliminated;  .  > 

-Since  y  =  x,  then  x  -=  ydy/dx  and  (3)  can,  be -written  as- 

^  F(y)  --x 

dx  - y  ■  ■  ■ 

F(y)  =  y  -  y^/3. 


(4) 
where 

(5) 


Let  us  show  that  any  trajectory  of  (4)  tends  to  a  periodic  orbit.  In 


Fig.  2  we  have  drawn  the  "fundamental”  curve  x 


starting  at  some  point...? 


FIGURE  12: 


Nb’te*  ^'fi'rs  t'^thht 

r-T.S'-:'  v 

and  .suchVithht 

'■  .  S' 

-thatv.s^y  \,<v'  0 . 


’.at”,  aIi>;?;.poiht^.^q;i,the6?l‘dft';;0'f''’' ‘the; 'curve 


.  y  O^'iand'-'  at  .^all  „■  po'ints-v 1 6  -".the't'.'right .;  o  f  ■• .  -:-x  -  '=,' 

..susvC 

.Ve  :have  iridicaj^ed  the,  '^arjt^s  of...the’-pl.ane  where 
by  "the  signs";4'*-+  '  and.  p-p  .’':re,spectively.' 


.?.( y)  ■■'•’■  ^n'd  sue  h 
■.V'dy/dX'..>,  0 


Note  also  that  'since 


■dx  yd.tj''f- ithe'  value’"'  of  ^.x  ,  ■?’;dn;5'ariy;,.'trajectory';,'raust;''lncre'ase  ih  the'  'upper 
half  ■'plane ‘’’an.dj,;deCrease.‘;in-, the  l'oWrer!'hal'f  ■■plane'jv  ''We'  have  .'.-indicated  •  this- 
'fact  'by- .the  .horizontal  'arrows '■!'  ,.'t... -tt I'--...  ■.'■'■ '  ''''v:!';'  ''  ' 


Consider’  -.a  ^traj e'ctdryljstartih'g'  ati'a  ’ p6Tht.'.r’  P ’  .r.  a.t  'some  distance  ■  from 


From ’-,(-'4)’ ’  because;  -e;  '.is'’ small,  the  value  of  dy/dx 


imst  be  large  and  negative.  Thus  the  trajectory  must  dr9p.very  rapidly 
tdwards’-the  fundament^.! . curve  x  =  F(y)’. .  The  trajectory  cannot  cross-  • 

.  '  •  the  •  fundamental  curve  for',  if  it  did,  the '  trajectory  would  be  in  a  region 

where  its’ slope  is  positive,  the  value  of  x  is  increasing  (watch  the 
■  arrow! )  and  the  value  of  y  is  decreasing.  The  trajectory  continues 
close  to  the  fundamental  curve  -until  it  passes  the  maximiun  point  B  and 
then,  as  can  be  seen  from  (4) ,  dy/dx  becomes  large  again  and  the 
trajectory  drops  rapidly  to  point  C.  At  C  ■  the  -trajectory  crosses  th’e 
fundamental  curve  and  dy/dx  becomes  negative .  The  trajectory  must’- stay 
close  to- the  fundamental  curve  for,  otherwise,  .dy/dx  becomes  very  large 
-and  the  trajectory  is  forced  back.  Note  also  that  the  trajec-tory  cannot 
.  ■  cross-'the  fundamental  curve  between  C  and  .  D  'because  above  CD’  the  • 

•  ■(  >  trajectory  must  ;go  toward  the  x-axis-with  a  posi-bive  slope  and  with  the 

value  of  X  decreasing,  conditions -which'  are  impossible  to  satisfy.  The 
'■■..■.conclusion  is.  then  'tha.-b- ’the  trajectory  is  clpse  to.  the  -fundamental  curve 
-  along  .CD  and ’then  rises  very  quickly  to  some  .’point.  E-  on  the  fundamental 
■1’,  curve.  ^  After  E  the  trajectbiy .keeps. -going  around  in  essentially  the  .-. 

'  ’_  -  -sa-rn'e  Orbit’-.  -  A 'rigorous .'proof.-that- a  periodic  trajectory  is  approached 
-’■  -  ■  -.'.can' be.  found’ -in  Stokes',  -J'Noniinear  -’Vibrations"  ,’ Appendix  IV. 

_■  .  ‘  ■  -  The  above  ’di’scussion  has  -  shown  that  .- 

-  the  .periodic  trajectory' can  be’.appro’xi-  .  - 

mated  by  two  .  arcs.  AB  and  CD  of  the.  _  '  • 
charac-teristlc  curve  and  two- -vertical 

•  •  1’ 

lines  BC  and  ’  DA  as  indicated  in  Fig.  3- 
To  this  approximation  we. may  find  the 

o 

period  T.  On  AB  and  DC  we  have 


FIGURE  3 
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X  =  F(y)  ■  but  on  BC  and  AD  we  have  x  _  constant.'.  Since  dt'  =  dj^^y. 


we  have 


T  =  jf=2j®  f  =  2(lny--V)'^ 


=  .  3 .  -  2  ln2  . 


This  naive  attempt  surprisingly  gives  the  correct  first  order  term 
for  the  period.  We  shall  see,  however,  that  to  get  a  correction  term  . 
requires  a  good  deal  of  analytic  manipula'tiqn. .  Let  us  suppose^"  that,  as 
indicated  in  Fig.  4,  we  have  a  periodic  trajectory  which  begins  at  some 
point  A^  on  the  fundamental  curve  and  goes  through  the  half  cycle 
indicated  by  the  points  . . . ,Arj.  The  coordinates  of  A^  will  be 

denoted  by  (x^jy^)  ,  1  <  i  <  V.  We  specify  x^  =  2/3,  thus  A2  is 


vertically  above  B  and  y^  =  1,  thus  A^  i  s  on  the-  same  horizontal 
line  as  B.  The  other  points  will  "be  specified  later. 

*  The  method  used  here  follows  very  closely  a  paper  by  W.  Wasow  in  . 
"Procedures  of  Syraposiiim  on  Non-Linear  Circuit  Theory,"  p.  87-97. 
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•  ■r*t  ’  ’  •»*', 


.then,- using  (4)  and  •  (5)-,  ;we'find.,^45%'^  -’I---.  '  •;*  v.’  -Vi:  •  v'  ..j.-.' 


' r  -  aSBCT  •■  -.k-"-  K.4  W-.  ,f' 

v4vv.  ■'  ■;■ 


On  the  arc  A- ' S^ce  ^^0  and  du/dy  <  0 

at  A,  and  since.  albng:,t  A^A^Ai '^'’^s^ha'^  <  0,  then  du  >  0  near  .  .  ■. 

■.  -  --A.  ^  :> 


,:•*••  fiA  ^  f 


If  du  stays  :;pDs.itive  in  going  from  A^  to.  A^,^,  .then  we  liras t  '•"  ’  _ 
^ave  in(7)'  "  '  '  i  '  '  ..  .... 


u  ...< 


’  'V---. . . 


If  du  =  0  at  Soirie  ' point..  .tnen"-  -uCyj-^i-^^^^'has^a.jmdximuiri^ there  and'‘ve  ”•' 

.  •;.•  .  ..•••■  •'  ••.^<v‘.  'v;.-vA-  •  '  ^.■‘;  •.■•»-■  •  . '€> 

have  for  y  <  y-^  ,..-t  :  7fv-i50  V\<«4  "  .'■/ 


'^(y)  - < 'u(y*)  -  |y-., ....  .  <'-''A-£_-x.  j...  ■  ,/y.-..; 

.  ■'■  ‘  ■  ■•■■• ....  -y^.  -  .‘1  .I-"  .i.y^^--..l:sv> .wy .  r>i^i...:  -  -  ■  'fw; ..-.  .  ‘ 

.....  .  ■•;..»«.....  .'t-V.  rs^iJvV'...'. .  .-“v'V'tv#  't  -  .  ■:  . 

■  ■  '  '•  ••».  -4.'  -  .*•  '•••  V  v’  psffc*  .  •  .  ■•*'*  *  *  '^!. 

.  Vs*.  6^.-.  '•  ‘4  ?•■’•' 3i  *‘’  ‘  * 

consequently  (8)  is  valid  for  allApoints. ' on .vi^a;A7^A ’  i'vp’. 

..•  •  .,  •  .•  ^  •-.'•A  W  '-f  t.  .»*N  •  > 

y."..  \  ' 

The  estimate  (8)  shows  that •  the' '.deviation. .fr.om,' the7.- arc  A^B  ,  is 
of.  the  order  e  as  long  as  the  ,traj.ectp'r.y'-'.is:j.h'ot  .’too  near..'B.’  To 
investigate  the  behavior  of  the  trajecto’i^  '.iii.  the  neighborhood  of  B, 


we  put 


y  F  1.  +  X  =  2/3  + 
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in  (4)..-  We  "find  "• 

(10) 


dT 


t  1  ,1/3,3  ^ 

1  + 


As  e  approaches  zero,  this  equation  may  be  approximated  by 
o 


(11)- 


dv  2 

—  +  V  +  T  =  0. 

dT 


.Equation  (ll)  is  a  Riccati  equation  and  with  the  help  of  the  standard 
transformation.'  . 


’  f- 


it  may  be  transformed  to  the  linear  second  order  equation 
(12) '  •  '  w”  +  TV  *=  0. 


The  solution, of  (12)  is 


®  .  . 


=  aA^(-T)  +  PB^(-t) 


where  -A^(t)  and  .•B^(t).  are  Airy  functions.  (For.  the  definitions  and 
properties  of  these  functions,  see  -Jeffreys,  "Methods. 'of  Mathematical 
.  Physics" . )  •  ■ 

The  appropriate  solution  of  (ll)  for  our  problem  is  the  one  that 

"matches"  with  the  solution  of  (4)  as  v  - >  oo  but  y  -  1  — >0.  Let 

us.  take  v  =  where  0  <  r  <  l/3.  Then  y  -  1  =  ^  — 5>  0  as 

e  — ^0.  From  (8)  and  (9)  we  find  that 

0<u=e'  (v  +T+e'  v/3).<e'  v 


I 
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This  ’simplifies  to  "  , 

/o  .  /  2  -  -1 

,  -  e  '  v/3  <  y  +  T  <  V 

“37  2  • 

and  for  v  =  -e  ,  we  see  that  v  + ’t  must  approach  zero  as  v  approaches 

'  ’  1/2 

infinity;  consequently,  v  is  approximately  equal  to  (-x)  '  for  large 
values  of  v.  Using  the  general  solution  (12) ,  we  have  that 


V  —  - 


aA^'(-x)  +  pB^'.(-x) 
■aA^(-T)  +  pB^(-x) 


and  by  studying  the  asymptotic  behavior  of  A^(-t)  and  B^(-t)  we  con- 

2  .  ■  • 

elude  that  if  v  +  t  =  0,  we  must  have 


V  -  - 


A^ ' ( -x) 
A^(-x) 


Formula  (13)  enables  us  to  complete  the  specification  of  the 
coordinates  of  k^k^k^.  The  point  A^  was  defined  by  =’2/3  or 
x^  =  0;  therefore, 


y„  =  1  -  e 


1/3 

Aj(0) 


The  point  A^  was  defined  by  y^  =  1  or  v  =  0.  If  x^  is  the  value 

of  T  such  that  A.'(-x_)  =  0,  then  ' 

1  3  ■ 

x^  =  2/3  ’+  . 


The  point  A,  has  not  been  defined  as  yet.  Let  us  take,  y .  =  1  -  e 
4-  •  ■  U 

where  0  <  r  <  l/3;  then  v  =  - >00  as  e  — >0.  Let  x.  be 

4 

the  smallest  value  of  'x  such  that  A.(-x.)  =  0.  (The. value  of  x,  is. 

14  -4 
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•  I 


approximately  2.33.)  Then  from  (13)  we  see  that  . 

(14)  •  "  =  2/3  + 

in  ordt r  to  ensure  that  v.  ’  approach  infinity.  '  .  •  '  • 

•  .  ®  /' 

The  point.  is  defined  as  the  intersection  of  the  trajectory  ■ 

with  thfe  x-axis;  consequently,  y^  =  O'.  To  find  ■  x^, we.  use  •• 

X,  -  x,  =  i  dx 

■  -  •  5  i  ..  -J.;  u  '■  . 

Trom  (4)-  We  may  express  u  as  follows: 

u  =  (y  -  1)^  ■+  J  (y  -  1)^‘  +  X  -  J  . 

Since  on  the  arc  the  value  of  x  is  always  increasing  (because 

dx  =  ydt  >  0.)  and  the  value  of  y  '  is  always  decreasing,  we  have 
X  -  2/3  >  0  and  |y- -  l|  <  1;  therefore, 

u  >  (y  -  1)-^[1  +  j(y  -  1)]  >  j(y  -  1)^  . 

Using  this  result,  we  find  that 


1-e  ,  „  ■  1-e 

/  3 


<  2  ^  ~  l)'’^dy  =.  0(e  "■);■ 

■  .0 


therefore , 


x^  -  x^  =  0(e^  ^)  =  oCe^'^^). 


The  point  .  is  defined  by  taking  y^^  =  y^  +  Again,  we  use  (4) 

for  dx  to  get  .  ,  .  ’ 


x^  -  x^  =  -  e 


r 

■’  ,  u 


y^+e 


V  I 

1 
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Noting  that  =  F(y^)  and  expressing  u  ln_  terms  of  coordinates  relative 

to  the  point  A^,  we  have 

^  .  -  •  • 

u  =  X  -  x^  +  |<y  -  y^)^  +  y^(y  -  y^)  ^  +  (y^  -  l)(y  -  y^) 

(15)  .  ....  .  • 

>  (y  -  y^)  [^y  -  y^)^  +  y^Cy  -  y^)  +  (y^  -  i)'] 

because  x  -  x^  is  positive  in  the  lower  half  plane.  The  bracket  in  (15) 
has  its  minimum  value  when 

y  -  y Y  =  -  3yY/2, 

that  is,  y  =  -y^/2.  'Since  this  value  of  y  Is  positive.,  we  conclude 
that  the  minimum  value  of  the  bracket  along  the  arc  4^4^  will  be  at  ^ 
y  =  0.  We  have  then  that 

u  >  (y  -  Vrj')  [j  y^  -  y^  y^’  1]  >  j  (y  -  y^)  • 

.  Using  this  estimate,  we  get 


(16) 
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(20) 


”t  =  -^  +  S'  ^ 

12  V  ■'  V  V 

■  yi  yi 


But 


(21) 

and 


^2  u  ^2 

f  ^  =  _2  +  C 

yi  ^  .^2  .y^  / 


*/  3 

""2  =  |-^^y2^  =  ^^^3^  -  ^  dF(y)'=  0.(e^/^). 


since  y^  -  y^  -  □(e^'^^)  .  We  have  .  • 


u 


:^  -  [1  +  a(e^/^)]  5  dF(y)  =  J  dF(y)  +’0(e) 

"^2  *  .  V 


(22) 


.y- 


J  ^(3-r-^)dF(y)  +  0(e)  =  ;^'^J^+6(e). 

’'2  ■  ■  ^2 


Also,  with  the  use'of  -(8),  we  see’ that 
(23) 


^2  y?  •  ••  ■.  ■  • 

^  ^  J  y~^  1  =  E  ’In-  (y  -  1) 

yi 


y^  y 


^2’ 


=  0( e  In  e) . 


Combining  (20) ,  (2l) ,  (22)  and  (23),  we  get 


(24) 


12 


f  Mz).+ 0(8^-^). 

y  '  ' 


Next,  we  evaluate 
T 


24 


(25) 


4 

i  ^  =  [1  +  0(e^)](x^  -■  x^) 


22/3. 


,2/3  + 


+  0(8'^/^  ^  ^)\ 
4 
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In  the  .same  way,  jus't  as  we  evaluated  x^-x^,  x^-x^  and  we  can 

ahow  that  .  ’  . 


"^47  .  ■  ^  .  u 

.  :  y. 


^  0(e^"^)  +‘0(e  In  e)  +  0(e) 


(26) 


Evaluating  (^A)  and  combining  with  (25)  and  (26),  we  get 

=  .8  +  3.5e^/^  +  o(e^/^) . 


Since  T  =  ^'^17’  obtain  finally 


T  =  1.6  +  7.0e^/^. 


